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Shackleton crater is nearly coincident with the Moon’s south pole. 
Its interior receives almost no direct sunlight and is a perennial 
cold trap 1-2 , making Shackleton a promising candidate location in 
[§] which to seek sequestered volatiles 3 . However, previous orbital and 
Earth-based radar mapping'' -8 and orbital optical imaging’ have 
yielded conflicting interpretations about the existence of volatiles. 
Here we present observations from the Lunar Orbiter laser 
Altimeter on board the Lunar Reconnaissance Orbiter, revealing 
Shackleton to be an ancient, unusually well-preserved simple crater 
whose interior walls are fresher than its floor and rim. Shackleton 
floor deposits are nearly the same age as the rim, suggesting that 
little floor deposition has occurred since the crater formed more 
l~3l than three billion years ago. At a wavelength of 1,064 nanometres, 
the floor of Shackleton is brighter than the surrounding terrain 
and the interiors of nearby craters, but not as bright as the interior 
walls. The combined observations are explicable primarily by 
downslope movement of regolith on the walls exposing fresher 
underlying material. The relatively brighter crater floor is most 
simply explained by decreased space weathering due to shadowing, 
but a one-micrometre-thick layer containing about 20 per cent 
surficial ice is an alternative possibility. 

Detailed study of the topography of Shackleton (Fig. la) offers the 
opportunity to improve understanding of processes that operate in 
permanently shadowed regions (Fig. lb). Crater geometry, age and 
preservation state are relevant for understanding the accumulation 
and preservation of volatiles as well as the processes that modify the 
lunar surface over geologic timescales. 

Our analysis uses observations from the Lunar Orbiter Laser 
Altimeter (LOLA) 10 , an instrument on NASA’s Lunar Reconnaissance 
Orbiter (LRO) mission. LOLA is a five-beam laser altimeter that operates 
at a wavelength of 1,064.4 nm with a 28-Hz pulse repetition rate. From 
LRO’s mapping orbit at ~50 km altitude, the instrument illuminates 
5-m-diameter spots on the lunar surface, returning up to 140 measure- 
ments of elevation per second; the five profiles enable characterization 
of bi-directional slopes over various baselines, and roughness from 
averaging of pulse elevations. In addition, from the spreading of 
backscattered laser pulses, LOLA obtains the root-mean-square 
(RMS) roughness of the surface within laser footprints. Finally, 
from the ratio of received to transmitted laser energy, LOLA measures 
the reflectance of the lunar surface at zero phase angle at the laser 
wavelength within laser spots. 

As of 1 December 2011, the LOLA instrument has accumulated 
more than 5.1 billion elevation measurements”. Because Shackleton 
lies nearly at a pole, where the LOLA coverage is densest, it is possible to 
construct a digital elevation model of unprecedented spatial resolution 
and radial accuracy. More than 5,000 LOLA tracks, referenced to the 
Moon’s centre of mass via precision orbits determined from radio 
tracking 12 aided by Earth-based laser tracking 13,1 ' 1 , were converted 


to topography. Track segments within the area of interest were geo- 
metrically corrected at orbit crossover points 15 . 

Figure la shows the topography of Shackleton crater sampled at 
10-m spatial resolution; individual measurements have an accuracy of 
~1 m with respect to the Moon’s centre of mass. The 40 km X 40 km 
topographic model of Shackleton is derived from 5.1 million elevation 
measurements with an average of 0.32 altimeter measurements in each 
10-m square area; the resolution is comparable to or better than other 
studies of Shackleton’s interior by images 9 . Earth-based radar 5 and 
orbital synthetic aperture radar 16 . The topography reveals the near- 
axisymmetric bowl-shaped nature of the crater, in which the crater rim 
and interior walls are well preserved. The depth/diameter ratio of the 
crater is 0.195 ± 0.025 (Table 1), which is consistent with other fresh 
simple craters 17 . 

Figure lc shows bi-directional surface slopes over 10-m baselines 
that quantify the uniform, steep inner walls; slopes approach the angle 
of repose, Slopes are greatest in the mid-levels of walls, which is in 
contrast to many crater walls on Mars, where near-vertically oriented 
cliffs of outwardly dipping coherent cap rock are exposed in the upper 
walls 18 . 

Surface roughness at a scale of 20-50 m is shown in Fig. Id. These 
data indicate that crater walls are smoother, within bounds of 
measurement uncertainty, at this spatial scale than the floor or rim, 
especially portions of lower walls aligned with mounds on the crater 
floor. Table 1 lists estimated roughness values of various crater com- 
ponents. Similarly, the average RMS roughness derived from spread- 
ing of individual laser footprints (Supplementary Information) is 
lower on the crater walls. Pulses returned from the steep walls of 
Shackleton are spread in time by >10 ns, but after correcting for the 
effect of local slope on a longer baseline, the pulse spreading due 
to surface roughness is somewhat less than on the crater floor or 
surrounding terrain. The floor can be divided into two regions, a flat 
portion and an elevated terrain. The roughness of tire mound unit 
increases at the largest scales due to its hummocky character, but it 
is smoother than the flat region at smaller scales, due to its paucity of 
craters. 

Figure 2 shows a more detailed view of the topography of 
Shackleton’s floor, which highlights the irregularly distributed deposits 
and numerous small craters (see also Fig. le). The largest mound of 
material has a relief of —210 m (Table 1) and the highest-local slope of 
any of the floor deposits is ~25°, which is below the angle of repose. 
Two areas of the floor show fan-shaped structures consisting of material 
that has been transported downslope from the interior walls in a 
manner commonly observed in craters of this size range 17 . The limited 
fan material around the margins of the crater floor, combined with the 
asymmetric distribution and slope properties of deposits, suggest 
that the predominant contribution to the fill is ejecta fallback with a 
secondary contribution from slumped wall deposits. 
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Figure 1 | Detailed characterization of Shackleton crater, a. Topography in 
km; b, percentage of time illuminated; c. 10-m baseline slopes in degrees; 
d, surface roughness shown as RMS residual in m; e, locations of crater counts 
used to determine relative ages; and f, zero-phase, 1,064-nm reflectance shown 
as II F. T opography, slopes and roughness are based on a 10-m spatial resolution 
grid of all available LOLA profiles. In a-d and f, x and y axes indicate spatial 
scale, where (0, 0) is the lunar south pole and colour scales show magnitude of 
plotted quantity. Panel e shows locations of craters counted to estimate relative 
age, plotted over 10-m slopes (colour coded as in inset). Crater regions in 
e correspond to: A, flat region of crater floor; A/B, entire crater floor; C, crater 


Table 1 | Parameters describing Shackleton crater 


Parjmeter 

Value 

Areocentric latitude of centre of rim (degrees) 

-89.655 

Areocentric longitude of centre of rim (degrees) 

129.174 

Lunar radius at floor centre (km) 

1,734.63 

Mean crater diameter at rim (km) 

21 

Mean depth, rim to floor (km) 

4.1 ±0.05 

Mean rim height above datum (km) 

1.3 

Range of floor topography (km) 
Area of crater at rim (krrr) 

-0.210 

-346 

Area of crater floor (km ? ) 

-38 

Estimated fill depth (km) 

-0.75 

Crater volume (km 2 3 ) 

640 ± 10 

Fill volume, including mounds (km 3 ) 

12 ± 1 

Maximum wall slope (degrees) 

35 

Average wall slope (degrees) 

30.5 

RMS roughness* of crater exterior (m) 

-1 

RMS roughness* of crater walls (m) 

<1 

RMS roughness* of crater floor (m) 

-1 

RMS roughness* of crater rim (m) 

-1 

l/F of crater exterior 

0.32 ± 0.04 

l/F of interior walls 

0.46 ±0.03 

l/F of interior floor 

0.43 ± 0.02 

Ratio of average depth/average rim diameter, d/0 

0.195 ±0.025 


See Fig. I legend for definition of l/F. 
• Within 5-m spot* 


wall; D, crater rim crest; E/F, inner rim annulus ( — 5.5 km); E, inner rim 
annulus excluding steep region (F); F, steep rim region within annulus; G, crater 
wall section; I, Shackleton crater deposits north of rim in flat areas; and X, 
secondary crater chains and clusters (removed from analysis). In f, reflectance 
is expressed as a radiance factor (J IF), which is defined as the ratio of the 
measured radiance / to the radiance F of an ideal diffusive surface in vacuum 
with 100% reflectance under the same illumination. Each dot represents a 
0.4 X 0.4 km pixel median average of LOLA’s spot 3 reflectance. Circles show 
21-km diameter of rim, 17-km diameter of the steepest portion of the walls, and 
7-krn diameter of the floor. 

Shackleton was previously assigned an Eratosthenian age 19 (middle 
lunar history; in the approximate interval 1 -3.2 Gyr before present) on 
the basis of its relatively fresh morphology, its lack of rays, and counts 
of superposed craters 3 using AMIE image data (50 m per pixel) and 
Arecibo radar data (20 m per pixel). Craters were counted within a 
crater diameter ( — 20 km) of the rim crest, avoiding obvious secondary 
craters. This analysis was subsequently revisited and resulted in an 
older, Imbrian age 3 (in the approximate interval 3.2-3.8 Gyr before 
present). 

Here we use a LOLA shaded relief map to advance previous work by 
individually dating different parts of the crater (see Fig. le and 
Supplementary Information) to investigate the processes that have 
operated since crater formation. LOLA observations permit dating 
of shadowed regions in the crater interior and allow spatially unbiased 
measurements of crater density due to uniformity in illumination 
conditions. In addition, illumination can be varied over the topographic 
model to enhance crater detection. On the basis of comparison of the 
several different areas of the rim of the crater, it is clear that the variable 
slopes of the rough crater rim have an influence on crater retention. For 
example, two areas of very flat terrain on the Shackleton flank within 
one crater diameter of the rim crest (I; Fig. 1 e) yield modelled crater ages 
of ~3.69 Gyr, whereas areas closer to the rim crest (within ~5.5 km) 
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Stereographic x axis (km) 

Figure 2 | High-resolution elevation map in stereographic projection of the floor of Shackleton. Elevations are contoured at 5-m intervals with colours 
indicating elevation with respect to 1,737.4 km. The axes indicate spatial scales. 


yield ages of ~1.21Gyr (F, steeper slopes; Fig. le) and ~2.91Gyr should be observed in terms of a deficit in the number of small craters. 

(E, fewer steep slopes; Fig. le). The flat areas of the Shackleton crater A rollover in the curve at small crater diameters is observed, but is very 

deposit (I;Fig.le)indicateanageof~3.69Gyr(SupplementaryTablel), similar to that seen outside the crater (I; Fig le and Supplementary 

older than the originally estimated age ofl. 3-3.3 Gyr(ref. 19)but close Fig. 1), and thus is more likely to be due to the typical destruction 

to the Upper Imbrian age of ~3.6Gyr estimated subsequently 3 , effects of superposed craters and other diffusive processes (such as 
Determination of ages of these individual regions permits quantitative micrometeorite bombardment or seismic shaking associated with 
investigation of how the crater has been modified. moonquakes triggered by stresses associated with impacts or tides) 

A critical question is the age of the permanently shadowed portions at these diameters. On the other hand, if volatiles were cold-trapped 

of the interior walls and floor of Shackleton. Analysis indicates that the by vapour diffusion into the regolith, as opposed to deposition in 

lunar spin axis has been at its approximate current orientation for surface layers, then muting of superposed craters may not have been 

~2 Gyr (ref. 20). If the crater had been accumulating volatiles in tire as significant. 

permanently shadowed areas over a period of this order, it is reasonable Figure 1 f shows profiles of 1 ,064-nm reflectance of Shackleton crater 

to hypothesize that its interior was resurfaced, covering and burying and its surroundings. A previous study 9 obtained images of the floor of 

craters and thus producing a younger relative age. Examination of the Shackleton from the Kaguya Terrain Camera at the time of maximum 

steep crater wall (C; Fig. le) yields a much younger crater retention age scattering illumination at the lunar south pole and observed no evid- 

of ~ 1 .44 Gyr, which could be consistent with either volatile mantling or ence for brightening; results were interpreted to indicate an absence of 

downslope mass wasting (Fig. lc) 21 . Examination of the permanently pure ice deposits on the crater floor. In the current study, LOLA pro- 

shadowed parts of the crater floor (A; Fig. le), the area where volatiles files assembled from numerous orbital passes at the most favourable 

plausibly accumulated, reveals a crater retention age of ~3.60 Gyr, conditions for obtaining reliable measurements of reflectance show 

essentially identical to the flat areas of the crater rim. Including the that the crater walls are anomalously bright relative to the surrounding 

rougher parts of the crater floor (B; Fig. le) produces an age of terrain. Similarly to other impact craters in this size range, this bright- 

~3.29 Gyr, which almost certainly reflects the influence of downslope ness could be due to downslope movement of material caused by 

transport, as observed on the inner part of the Shackleton rim, on the micrometeorite and small projectile bombardment on steep slopes, 

retention of craters. The similar age of the Shackleton crater exterior or by seismic shaking. The cascading of regolith material downslope 

and floor (~3.69 and 3.60 Gyr) is evidence that if volatiles accumulated exposes optically less mature surfaces than those developed and 

in this cold trap for ~ 10 9 years, they were not in sufficient quantity to retained on lower slopes. 

alter in a statistically significant sense the size- frequency distribution of At a wavelength of 1 ,064 nm, the floor of Shackleton crater is darker 
superposed craters in the size range counted. than its interior walls, but both floors and walls are considerably 

Knowledge of the crater size-frequency distribution (Supplemen- brighter than the surrounding terrain, including the interiors of nearby 

tary Fig. 1) permits an assessment of the minimum amount of depos- craters that are both shadowed and sunlit. The relative brightness of the 

ition that could occur without disrupting this distribution. The craters floor relative to surroundings requires explanation. Micrometeorite 

counted on Shackleton’s floor ranged from ~250 m up to ~500 m in bombardment and impingement of the solar wind produce ‘space 

diameter, representing a range of fresh crater depths of ~50-100m weathering’ of exposed geologic materials that reddens and darkens 

(ref. 22). Accumulation of volatiles to thicknesses in excess of 20-50 m their surfaces 23 . The former (bombardment) would be less significant 

would significantly alter the crater size-frequency distribution and in the shadowed interior of Shackleton than in the surrounding region 
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because the interior of the crater has not been exposed to the Sun for 
more than 2 Gyr (ref. 20), whereas the latter (solar wind impingement) 
might be enhanced in permanently shadowed craters 3 ' 1 . Thus the floor 
brightness enhancement could be explained by a dearth of space 
weathering by micrometeorite bombardment. 

Volatile deposition is an alternative possibility. Under the conser- 
vative assumption that water ice has a 1,064-nm reflectance twice that 
of the lunar regolith 35 and that both are observed at zero phase, the 
measured reflectance of the floor can be explained by a micrometre- 
thick surface layer (the depth over which the laser backscatter 
measurement is sensitive) of 22% ice mixed with rock 36 . Greater ice 
contents distributed throughout a thicker layer are possible but cannot 
be constrained from LOLA’s reflectivity measurement. For com- 
parison, far-ultraviolet reflectance of permanently shadowed regions 
from the LRO Lyman Alpha Mapping Project (LAMP) is consistent 
with ~l-2% surface water frost 37 . 

Results from the LRO Mini-RF orbital radar for the interior of 
Shackleton 16 provide additional insight. Regions with thick ice deposits 
are expected to have circular polarization ratios (CPR) >1, but such 
high ratios can also be explained by surface roughness. High-resolution 
images from ground-based radar 6 show that some areas with high CPR 
lie inside the rim of Shackleton. Mini-RF data 16 reveal that CPR values 
decrease with depth within the crater, and CPR values on the floor of 
Shackleton crater are predominantly <1. The pixels with CPR values 
in excess of unity are distributed heterogeneously throughout the 
crater walls, correlate generally with regions of high roughness 
observed by LOLA (Fig. Id) and include some sunlit areas (Fig. lb). 
Although some contribution to high CPR values from volatiles is 
possible, and would imply a process in which volatile deposition 
operates very rapidly in comparison to the rate of removal, the com- 
bined data suggest that the higher floor reflectance is due primarily to 
the dearth of space weathering in this shadowed environment. 

In considering why Shackleton’s interior walls have a higher reflec- 
tance than its floor, it is instructive to note that wall brightening is not 
restricted to areas that are continuously shadowed but extends to the 
upper illuminated portions. Consequently, a higher concentration of 
surface volatiles on the walls than present on the floor is an unlikely 
explanation. More likely is downslope movement of regolith material 
on the steep crater walls that has exposed brighter underlying material; 
downslope movement is consistent with the observed slopes near the 
angle of repose and the roughness and morphology of Shackleton’s 
interior walls (Fig. lc, d), as well as with Earth-based radar backscatter 6 
and Mini-RF 16 measurements. 

METHODS SUMMARY 

LOLA, an instrument on board the LRO spacecraft, outputs five beams per laser 
pulse that are backscattered from the lunar surface and detected in the instru- 
ment’s receiver. The relevant measurement is the time of flight of each individual 
laser pulse, which can be converted to a range of the spacecraft to the lunar surface 
given knowledge of the position of the spacecraft with respect to the Moon’s centre 
of mass. Radial range errors were minimized by geometric adjustment of altimetric 
tracks within the study area. In practice, profiles of one-way range from the LOLA 
instrument to the lunar surface along the spacecraft ground track were converted 
to lunar radius at each bounce point using the reconstructed orbit of LRO. 
Topography was determined by subtracting a sphere of 1,737.4 km from each 
radius measurement. 

Slopes were calculated from a two-laser-spot fit at 20-50-m length scales, and 
RMS roughness represents a standard deviation about a plane fitted by least- 
squares to two laser shots along-track, from which at least four valid spots are 
returned out of a possible total of ten. The plane has dimensions of 90- 1 00 m in the 
longest axis and 10-40 m in the shortest axis, depending on the positions of the 
spots returned. 

The ratio of returned to transmitted pulse energy is a measure of surface reflec- 
tance at the laser wavelength. The transmitted and returned pulse energies were 
measured by integrating the area under the pulses. Observations of 1,064-nm 
reflectance were derived from LOLA tracks crossing Shackleton for days 130- 
149 in 2010. 


Full Methods and any associated references are available in the online version of 
the paper at www.nature.com/nature. 
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METHODS 

The Lunar Orbiter Laser Altimeter (LOLA), an instrument aboard the Lunar 
Reconnaissance Orbiter (I,RO) spacecraft, outputs five beams per laser pulse that 
are backscattered from the lunar surface and detected in the instrument's receiver. 
The relevant measurement is the time of flight of each individual laser pulse, which 
can be converted to a range of the spacecraft to the lunar surface given knowledge 
of the position of the spacecraft with respect to the Moon's centre of mass. 
The timing of LOLA instrument events was derived from the LRO ultrastable 
oscillator, which is monitored by ground tracking stations. Time systems on board 
LRO used Coordinated Universal Time (UTC) to correlate spacecraft Mission 
Elapsed Time (MET) to ground time. The analysis of LOLA data used 
Barycentric Dynamical Time as its primary time system. Spacecraft states relative 
to the Solar System Barycentre (SSB) at the laser transmit and detector receive 
times were projected along the instrument boresight and return path vectors to 
match the observed time of flight, correcting for the aberration of light and 
general-relativistic time delays. SSB states were determined in the Earth Mean 
Equator of 2000 ( J2000) inertial reference frame using lunar spacecraft trajectories 
and the DE421 planetary ephemeris 38 

During its lunar mapping mission, the LRO spacecraft is tracked using S-band 
Doppler and range data by the Universal Space Network, Deep Space Network and 
White Sands Missile Range. The precise reconstruction of LRO orbits used 
Doppler tracking observations from these stations as well as laser ranging to 
LRO 1 ' 1 from the Goddard Space Flight Center and participating members of the 
International Laser Ranging Service. Precision orbit determination was 
accomplished using the NASA/Goddard Space Flight Center's GEODYN system 
of programs 39 using the GLGM-3 gravity model 50 as a reference. GEODYN 
numerically integrates the spacecraft Cartesian state and force-model partial 
derivatives by employing a high-order Cowell predictor-corrector model. In addi- 
tion to a model of the lunar gravity field, the force modelling included point mass 
representations for the Sun and planets. Solar radiation pressure, measurement 
and timing biases, and tracking station coordinates were also estimated. 

Radial range errors were minimized by geometric adjustment of altimetric 
tracks within the study area. In practice, profiles of one-way range from the 
LOLA instrument to the lunar surface along the spacecraft ground track were 
converted to lunar radius at each bounce point using the reconstructed orbit of 
LRO. Topography (Fig. la) was determined by subtracting a sphere of 1,737.4 km 
from each radius measurement. Topographic measurements were binned and 
interpolated within a 40 km X 40 km area at 10-m spatial resolution. 


The LOLA digital elevation model (DEM) used in this analysis was combined 
with a lunar ephemeris 51 to characterize the solar illumination conditions of 
Shackleton and surroundings (Fig. lb). A polar gnomonic projection on which 
great circle paths plot as straight lines was applied to the DEM to calculate lighting 
conditions that are accurate over geological timescales. 

The confluence of orbit ground tracks in the vicinity of the lunar poles, com- 
bined with LOLA's multi-beam profiling capability, permitted slopes over a range 
of baselines and directions to be determined. In this study, slopes (Fig. 1c) were 
calculated from a two-laser-spot fit at 20-50-m length scales. In similar fashion, 
RMS roughness (Fig. Id) was calculated as a standard deviation about a plane 
fitted by least-squares to two laser shots along- track, from which at least four valid 
spots were returned out of a possible total of ten. The plane had dimensions of 
90-100 m in the longest axis and 10-40 m in the shortest axis, depending on the 
positions of the spots returned. An independent measure of surface roughness, 
discussed in Supplementary Information and shown in Supplementary Fig. 2, used 
the spreading in time of backscattered pulses. The spreading of LOLA’s 
backscattered pulses provides a measure of the RMS roughness of the surface at 
a smaller scale— the 5-m diameter of the laser footprints on the lunar surface. 
As with plane-deviation roughness shown in Fig. Id, the pulse-spread-derived 
roughness featured a correction for local slopes. 

The ratio of returned to transmitted pulse energy is a measure of surface 
reflectance at the laser wavelength of 1,064 nm. The transmitted and returned 
pulse energies were measured by integrating the area under the pulses. LOLA’s 
measurement of reflectance is calibrated only in a relative sense, with respect to 
pre-launch testing, as the instrument lacks a source with known brightness in 
flight. Observations of reflectance were derived from LOLA tracks crossing 
Shackleton and environs from day of year 130-149 2010, which represented the 
most favourable time period for stable reflectance measurements due to the 
geometry of spacecraft terminator crossings. 
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